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2. COMMUNICATION: 



Due to a clerical error, the international application has not been published promptly after the 
expiration of 18 months from the priority date, as provided in Article 21(2) of the PCT. 

Consequently, the international publication will only take place on 01 November 2001 
(01 .1 1.01). Meanwhile, the International Bureau (WO) will communicate the international 
application to each designated Office, in accordance with Article 20. 
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Applicant's or agent* s file reference 
P23847 A/ER A/P P P 


See Notification of Transmittal of International 
FOR FURTHER ACTION Preliminary Examination Report (Form PCT/IPEA/416) 


International application No. 
PCT/GB00/01675 


International filing date (day/month/year) 
02/05/2000 


Priority date (day/month/year) 
01/05/1999 


International Patent Classification (IPC) or national classification and IPC 
G06F17/00 


Applicant 

THE COURT OF NAPIER UNIVERSITY et al. 





1 . This international preliminary examination report has been prepared by this International Preliminary Examining Authority 
and is transmitted to the applicant according to Article 36. 

2. This REPORT consists of a total of 12 sheets, including this cover sheet. 

□ This report is also accompanied by ANNEXES, i.e. sheets of the description, claims and/or drawings which have 
been amended and are the basis for this report and/or sheets containing rectifications made before this Authority 
(see Rule 70.16 and Section 607 of the Administrative Instructions under the PCT). 

These annexes consist of a total of sheets. 



3. This report contains indications relating to the following items: 



Basis of the report 



Reasoned statement under Article 35(2) with regard to novelty, inventive step or industrial applicability; 
citations and explanations suporting such statement 

Certain documents cited 

Certain defects in the international application 

Certain observations on the international application 



II 


□ 


III 




IV 


□ 


V 


IS 


VI 


□ 


VII 




VIII 





Date of submission of the demand 
30/11/2000 


Date of completion of this report 
17.08.2001 


Name and mailing address of the international 
preliminary examining authority: 

^ European Patent Office 
/§)) D-80298 Munich 

C^' Tel. +49 89 2399 - 0 Tx: 523656 epmu d 
Fax: +49 89 2399 - 4465 


Authorized officer 

Barba.M f ^ j) 

Telephone No. +49 89 2399 2732 X^g*-*^ 



Form PCT/IPEA/409 (cover sheet) (January 1994) 




INTERNATIONAL PRELIMINARY 

EXAMINATION REPORT International application No. PCT/G BOO/0 1675 



I. Basis of the report 

1 . With regard to the elements of the international application (Replacement sheets which have been furnished to 
the receiving Office in response to an invitation under Article 14 are referred to in this report as "originally filed" 
and are not annexed to this report since they do not contain amendments (Rules 70. 16 and 70. 17)): 
Description, pages: 

1 -40 as originally filed 

Claims, No.: 

1 -40 as originally filed 

Drawings, sheets: 

1/14-14/14 as originally filed 



2. With regard to the language, all the elements marked above were available or furnished to this Authority in the 
language in which the international application was filed, unless otherwise indicated under this item. 

These elements were available or furnished to this Authority in the following language: , which is: 

□ the language of a translation furnished for the purposes of the international search (under Rule 23.1 (b)). 

□ the language of publication of the international application (under Rule 48.3(b)). 

□ the language of a translation furnished for the purposes of international preliminary examination (under Rule 
55.2 and/or 55.3). 

3. With regard to any nucleotide and/or amino acid sequence disclosed in the international application, the 
international preliminary examination was carried out on the basis of the sequence listing: 

□ contained in the international application in written form. 

□ filed together with the international application in computer readable form. 

□ furnished subsequently to this Authority in written form. 

□ furnished subsequently to this Authority in computer readable form. 

□ The statement that the subsequently furnished written sequence listing does not go beyond the disclosure in 
the international application as filed has been furnished. 

□ The statement that the information recorded in computer readable form is identical to the written sequence 
listing has been furnished. 

4. The amendments have resulted in the cancellation of: 

□ the description, pages: 

□ the claims, Nos.: 
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□ 



the drawings, 



sheets: 



5. □ This report has been established as if (some of) the amendments had not been made, since they have been 

considered to go beyond the disclosure as filed (Rule 70.2(c)): 

(Any replacement sheet containing such amendments must be referred to under item 1 and annexed to this 
report.) 

6. Additional observations, if necessary: 

III. Non-establishment of opinion with regard to novelty, inventive step and industrial applicability 

1. The questions whether the claimed invention appears to be novel, to involve an inventive step (to be non- 
obvious), or to be industrially applicable have not been examined in respect of: 

□ the entire international application. 
E*3 claims Nos. 16-32, 38-40. 



□ the said international application, or the said claims Nos. relate to the following subject matter which does 
not require an international preliminary examination {specify): 

the description, claims or drawings (indicate particular elements beloW) or said claims Nos. 1 6, 38 are so 
unclear that no meaningful opinion could be formed (specify): 
see separate sheet 

□ the claims, or said claims Nos. are so inadequately supported by the description that no meaningful opinion 
could be formed. 

Kl no international search report has been established for the said claims Nos. 17-32, 39-40. 

2. A meaningful international preliminary examination cannot be carried out due to the failure of the nucleotide 
and/or amino acid sequence listing to comply with the standard provided for in Annex C of the Administrative 
Instructions: 

□ the written form has not been furnished or does not comply with the standard. 

□ the computer readable form has not been furnished or does not comply with the standard. 

V. Reasoned statement under Article 35(2) with regard to novelty, inventive step or industrial applicability; 
citations and explanations supporting such statement 

1. Statement 

Novelty (N) Yes: Claims 3-15,35-37 
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No: 



Claims 



1 , 2, 33-34 



Inventive step (IS) 



Yes: 
No: 



Claims 
Claims 



3-15, 35-37 



Industrial applicability (IA) 



Yes: 
No: 



Claims 
Claims 



1-15, 33-37 



2. Citations and explanations 
see separate sheet 

VII. Certain defects in the international application 

The following defects in the form or contents of the international application have been noted: 
see separate sheet 

VIII. Certain observations on the international application 

The following observations on the clarity of the claims, description, and drawings or on the question whether the 
claims are fully supported by the description, are made: 
see separate sheet 
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Reference is made to the following documents: 

D1 : WO 96 08992 A (SHOSHAN HERBERT Z ;UNIV RAMOT (IL); AKSELROD 
SOLANGE (IL); KESELBR) 28 March 1996 (1996-03-28) 

D3: SAVA H ET AL: 'APPLICATION OF THE MATCHING PURSUIT METHOD 
FOR STRUCTURAL DECOMPOSITION AND AVERAGING OF 
PHONOCARDIOGRAPHY SIGNALS' MEDICAL AND BIOLOGICAL 
ENGINEERING AND COMPUTING, GB,PETER PEREGRINUS LTD. 
STEVENAGE, vol. 36, no. 3, 1 May 1998 (1998-05-01), pages 302-308, 
XP000751653 ISSN: 0140-0118 

The following document (D) was not cited in the international search report. A 
copy of the document is annexed to the Written Opinion and the numbering will be 
adhered to in the rest of the procedure: 

D2: Proceedings of Computers in Cardiology, IEEE, September 23-26 1991 , 
Venice Italy, pages 393-396, D. Morlet et al: "Time-Scale Analysis of High 
Resolution Signal Averaged Surface ECG Using Wavelet Transformation" 



0 With regard to present claims 17 to 32 and 39 to 40, it is noted that no 
international search report has been established in respect of the above 
mentioned set of claims. 

Consequently this International Preliminary Examining Authorithy, under the 
provisions of Rule 66.1 (e), does not need to carry out an international preliminary 
examination in respect to the subject matter of present claims 17 to 32 and 39 to 
40. 

0.1 This International Preliminary Examination will be therefore limited to the subject 
matter of present claims 1 to 16 and 33 to 38. 
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Re Item III 

Non-establishment of opinion with regard to novelty, inventive step and 
industrial applicability 

1 When considering the extent of the clarity problems regarding present 

independent claims 16 and 38 (see Item VIII below) this International Preliminary 
Examining Authority considers as not possible to give an opinion as to novelty, 
inventive step and industrial applicability in respect of the above mentioned claim 
16 and 38. 



Re Item V 

Reasoned statement under Rule 66.2(a)(ii) with regard to novelty, inventive 
step or industrial applicability; citations and explanations supporting such 
statement 

2 Insofar as the present text of independent claim 1 can be understood in the light 
of the description (see Item VIII below), it appears that the subject matter of claim 
1 does not fulfill the requirements of novelty as set out in Article 33 (2) PCT, the 
reasons therefor being the following. 

2.1 Document D1 , that provisionally is considered as representing the closest prior 
art, discloses (see from page 9 line 27 to page 12 line 12; from page 13 line 29 to 
page 15 line 12) a method to analyse an ECG signal by using wavelet 
decomposition of said ECG signal. 

Thus, the subject matter of claim 1 is not novel against the method known from D1 
(Article 33 (2) PCT). 



3 Dependent claim 2 does not contain any features which, in combination with the 
features of any claim to which it refers, meet the requirements of the PCT in 
respect of novelty, because the feature of using discrete wavelet transforms is 
also disclosed in document D1 (Article 33 (2) PCT). 
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4 Insofar as the present text of dependent claims 3 to 15 can be understood in the 
light of the description (see Item VIII below), it also appears that the subject matter 
of claims 3 to 15 does not involve an inventive step and therefore the 
requirements of Article 33 (3) PCT are not met, the reasons therefor being the 
following. 

4.1 Document D2 discloses (see page 393 left column lines 1 to 21 ; from page 393 
right column line 16 to page 394 left column line 1 ; page 394 right column lines 1 
to 23; figure 4 and figure 5) a method to analyse an ECG signal using continuous 
wavelet decomposition including the following features: 

i) computing wavelet energy surfaces of said ECG signal and plotting said 
wavelet energy surfaces against parameters a and b of the wavelets bases; 

ii) constructing a contour plot and a surface plot to display said wavelet 
decomposition of said ECG signal; 

iii) constructing a 2D or 3D energy scalogram to display said wavelet 
decomposition of said ECG signal. 

4.1a Moreover, the subject matter of dependent claims 7, 8, 9, 10, 1 1 and 15 are 
considered obvious. 

5 Therefore, the subject matter of dependent claims 3 to 15 does not include an 
inventive step contribution in respect of the method known from the combination 
of the method of D1 and the method of D2 (Article 33 (3) PCT). 

6 The below mentioned lack of clarity notwithstanding (see Item VIII below), the 
subject-matter of claim 33 is not novel in the sense of Article 33 (2) PCT because 
document D3 discloses (see from page 303 left column line 1 1 to page 305 right 
column line 54) a method to decompose a cardiac signal using a matching pursuit 
algorithm. 

7 The below mentioned lack of clarity notwithstanding (see Item VIII below), the 
subject-matter of independent apparatus claim 34 is not novel in the sense of 
Article 33 (2) PCT for the same reasons, mutatis mutandis already mentioned in 
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paragraph 2.1 of this International Preliminary Examination Report. 

8 The below mentioned lack of clarity notwithstanding (see Item VI 1 1 below), the 
subject-matter of dependent claims 35 to 37 is not inventive in the sense of Article 
33 (3) PCT for the same reasons, mutatis mutandis, already mentioned in 
paragraph 3 to 4.1b of this International Preliminary Examination Report. 

9 With regard to the assessment of the present claims 1 to 15 and 33 to 37 on the 
question whether they are industrially applicable, the following is stated. 

The below mentioned lack of clarity notwithstanding (see Item VIII below), it 
appears that the subject matter of present claims 1 to 16 and 33 to 37 relates to a 
method and apparatus to decompose waveforms of a cardiac signal, therefore it 
fulfills the requirements of industrial applicability as set out in Article 33 (4) PCT. 



Re Item VII 

Certain defects in the international application 

10 At page 6 line 28 of the description the wording "me" should be amended as 
"method". 

At page 9 line 17 of the description the wording "scrologram" should be amended 
as "scalogram". 

10.1 Present independent claims are not in the two-part form in accordance with Rule 
6.3(b) PCT, which in the present case would be appropriate, with those features 
known in combination from the prior art (document D1) being placed in a 
preamble (Rule 6.3(b)(i) PCT) and with the remaining features being included in a 
characterising part (Rule 6.3(b)(ii) PCT). 
Independent claims should therefore be redrafted accordingly. 
In addition, it should be clear from the description which features of the claimed 
subject-matter are known from documents D1 , D2 and D3 (see the PCT 
Guidelines PCT/GL/3 III, 2.3a). 
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10.2 Contrary to the requirements of Rule 5.1 (a)(ii) PCT, the relevant background art 
disclosed in the documents D1, D2 and D3 is not mentioned in the description, nor 
are these documents identified therein. 

10.3 The features of the claims are not provided with reference signs placed in 
parentheses (Rule 6.2(b) PCT). 

10.4 Furthermore, at page 40, last paragraph, the description contains general 
statements that the extent of protection may be expanded in some vague and not 
precisely defined way. Such general statements shall be deleted as contrary to 
Article 6 PCT, cf. also PCT Preliminary Examination Guidelines, C-lll, 4.3a. 



Re Item VIII 

Certain observations on the international application 

1 1 Present independent claim 1 is not clear and as such it does not fulfill the 
requirements of Article 6 PCT for the following reasons. 

11.1 The wording of present claim 1 is too broad and vague and it does not enable the 
person skilled in the art to carry out the invention without any further inventive 
effort, which is against Article 6 PCT. 

It is the opinion of this Authority that the person skilled in the art would be unable, 
on the basis of the information given in the application as filed, to extend the 
particular teaching of the description to the whole of the fields claimed by using 
routine methods of experimentation or analysis. Therefore, the wording of present 
claim 1 should be amended in order to properly limit the extent of the subject 
matter claimed in accordance with the subject matter as disclosed in the 
application as a whole. 

12 Dependent claim 2 is also vague and unclear and therefore it does not meet the 
requirements of Article 6 PCT, the reasons therefor being the following. 
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12.1 The wording "comprising the step of employing discretized wavelet transform 
analysis to process the cardiac waveform" used in claim 2 is vague and unclear 
and leaves the reader in doubt as to the meaning of the technical features to 
which it refers, thereby rendering the definition of the subject matter of said claim 
unclear (Article 6 PCT). 

13 Dependent claim 3 is also vague and unclear and therefore it does not meet the 
requirements of Article 6 PCT, for the same reason already mentioned in 
paragraph 12.1 of this International Preliminary Examination Report and for the 
following additional reason. 

13.1 The wording "discretized continuous" used in claim 3 does not have a clear 

technical meaning and as such it leaves the reader in doubt as to the meaning of 
the technical features to which it refers, thereby rendering the definition of the 
subject matter of said claim unclear (Article 6 PCT). 
This is also so because a wavelet basis either is discrete or is continuous; 
moreover it appears from the description that the method of the application is 
specified in base of continuous Morlet wavelet basis functions. 

14 Dependent claim 7 is also vague and unclear and therefore it does not meet the 
requirements of Article 6 PCT, the reasons therefor are the following. 

14.1 The wording "derive the cardiac signal" used in claim 7 is unclear and as such it 
leaves the reader in doubt as to the meaning of the technical features to which it 
refers, thereby rendering the definition of the subject matter of said claim unclear 
for the following reasons. 

The wording of present claims 1 to 6 does not specify whether the named "cardiac 
signal" is an analogical signal or it is a digitized signal, while from the wording of 
dependent claim 7 it appears that what the applicant meant was a digitized 
cardiac signal. This fact creates in the reader a state of uncertainty as to the 
extent of the subject matter claimed, which is against the provisions of Article 6 
Pet. 
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15 Dependent claims 10 and 1 1 are also vague and unclear and as such they do not 
fulfill the provisions of clarity of Article 6 PCT for the following reasons. 

15.1 The wordings "coherent structures" and "for clinical use" used in claim 10 and 1 1 
respectively do not have a clear technical meaning and as such they leave the 
reader in doubt as to the meaning of the technical features to which it refers, 
thereby rendering the definition of the subject matter of said claim unclear (Article 
6 PCT). 

16 Claim 15 does not meet the requirements of Article 6 PCT in that the matter 
claimed is not clearly defined. The claim attempts to define the subject matter in 
terms of the result to be achieved which merely amounts to a statement of the 
underlying problem. The technical features necessary for achieving this result 
should be added. 

17 Independent claim 16 does not meet the requirements of Article 6 PCT in that the 
matter claimed is not clearly defined. The overall wording of claim 16 does not 
provide the required intelligibility in order to be understood by the average skilled 
person. It is necessary to redraft the claim using a more understandable manner, 
in order to meet the requirements of Article 6 PCT (clarity in the sense of 
intelligibility). 

18 Present independent claim 33 is not clear and as such it does not fulfill the 
requirements of Article 6 PCT for the same reasons already mentioned in 
paragraph 11.1 of this International Preliminary Examination Report. 

19 Independent claim 34 and dependent claims 35 to 37 are also vague and unclear 
and as such they do not fulfill the requirements of Article 6 PCT for the following 
reasons. 

19.1 Claims 34 to 37 do not meet the requirements of Article 6 PCT in that the matter 
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claimed is not clearly defined. The functional statements included in the wording 
of present claims 34 to 37 do not enable the skilled person to determine which 
technical features are necessary to perform the stated functions. 




20 Independent claim 38 does not meet the requirements of Article 6 PCT in that the 
matter claimed is not clearly defined for the same reasons as above mentioned in 
paragraph 17 of this International Preliminary Examination Report. 
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International Application No. 



International Filing Date 



02 



MAY 2000 
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Unittu Rangaom Patent Office ! 



Applicant's or agent 1 s file reference 

(if desired) (12 characters maximum) P23847A/JMK 



Box No. I TITLE OF INVENTION 

"Method of Analysis of Medical Signals" 


Box No. II APPLICANT 


Name and address: (Family name followed by given name; for a legal entity, full official 
designation. The address must include postal code and name of country. The country of the 
address indicated in this Box is the applicant s State (that is, country) of residence if no State 
of residence is indicated below.) 

The Court of Napier University 
Merchiston Tower 
Colinton Road 
EDINBURGH 
EH10 5DT 
United Kingdom 


| | This person is also inventor. 


Telephone No. 


Facsimile No. 


Teleprinter No. 


State (that is, country) of nationality: 

United Kingdom 


State (that is, country) of residence: 

United Kingdom 


This person is applicant ( 1 all designated 1 1 all designated States except 1 | the United States i | the States indicated in 
for the purposes of: 1 1 States |A | the United States of America 1 1 of America only 1 Ithe Supplemental Box 


Box No. Ill FURTHER APPLICANT^) AND/OR (FURTHER) INVENTOR(S) 


Name and address: (Family name followed by given name; for a legal entity, full official 
designation. The address must include postal code and name of country. The country of the 
address indicated in this Box is the applicant s State (that is, country) of residence if no State 
of residence is indicated below.) 

ADDISON, Paul Stanley 
58 Buckstone Crook 
EDINBURGH 
EH10 6UR 
United Kingdom 


This person is: 

| | applicant only 

|y | applicant and inventor 

1 inventor only (If this check-box 
is marked, do not Jill in below.) 


State (that is, country) of nationality: 


State (that is, country) of residence: - 

United-Kingdom 


United Kingdom 




This person is applicant | | all designated I | all designated States except the United States I I the States indicated in 
for the purposes of: I I States I 1 the United States of America LcJ of America only 1 1 the Supplemental Box 


|Xl Further applicants and/or (further) inventors are indicated on a continuation sheet. 


Box No. IV AGENT OR COMMON REPRESENTATIVE; OR ADDRESS FOR CORRESPONDENCE 


The person identified below is hereby/has been appointed to act on behalf rrn i i representative 
of the applicants) before the competent International Authorities as: LCJ ° 1 1 r 


Name and address: (Family namefollowed by given name; for a legal entity, full official 
designation. The address must include postal code and name of country.) 

Murgitroyd & Company 
373 Scotland Street 
Glasgow 
G5 8QA 
United Kingdom 


Telephone No. 

0141 307 8400 


Facsimile No. 

0141 307 8401 


Teleprinter No. 


1 | Address for correspondence: Mark this check-box where no agent or common representative is/has been appointed and the 
1 1 space above is used instead to indicate a special address to which correspondence should be sent. 
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Sheet No. 



Continuation of Box No. HI FURTHER APPLICANT(S) AND/OR (FURTHER) INVENTOR(S) 



If none of the following sub-boxes is used, this sheet should not be included in the request 



Name and address: (Family name followed by given name; for a legal entity, full official 
designation. The address must include postal code and name of country. The country of the 
address indicated in this Box is the applicant s State (that is, country) of residence if no State 
of residence is indicated below.) 

WATSON, James Nicholas 
34 Fowler Terrace 
EDINBURGH 
EH11 1DA 
United Kingdom 



This person is: 

| | applicant only 

\)C\ applicant and inventor 

| [ inventor only (If this check-box 
is marked, do not Jill in below.) 



State (that is. country) of nationality: 

United Kingdom 



State (that is, country) of residence: 

United Kingdom 



This person is applicant 
for the purposes of: 



j j all designated 



States 



□ all designated States except 
the United States of America 



fyi the United States 
I ^ I of America only 



□ the States indicated in 
the Supplemental Box 



Name and address: (Family name followed by given name; for a legal entity, full official 
designation. The address must include postafcode and name of country. The country of the 
address indicated in this Box is the applicant 's State (that is, country) of residence if no State 
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1 "Method of Analysis of Medical Signals" 
2 

3 This invention relates to a method of analysis of, 

4 medical signals, and in particular to a method of 
~~5 decomposition ot cardiac signals using wavelet 

6 transform analysis. Specifically the invention relates 

7 to an improved method of resuscitation of patients in 

8 cardiac arrest . 
9 

10 In the UK, coronary heart disease is the second 

11 greatest contributor to deaths of people under 75. The 

12 social and economic consequences of these death rates 
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2 

1 are enormous. The current survivability rates of 

2 patients after sudden cardiac failure are around 1:10. 
3 

4 Ventricular tachyarrhythmias, specifically ventricular 

5 fibrillation (VF) , are the primary arrhythmic events in 

6 cases of sudden cardiac death. Administration of 

7 prompt therapy to a patient presenting with such 

8 symptoms can however lead to their successful 

9 resuscitation. Until recently, the only indicators of 

10 likelihood of survival of a patient to hospital 

11 discharge were traditional variables such as emergency 

12 service response time or bystander cardio -pulmonary 

13 resuscitation (CPR) . 
14 

15 In most cardiac complaints, analysis of a surface 

16 electrocardiogram (EKG). of the presenting patient is a 

17 rich source of information. However, until recently, a 

18 surface EKG recorded during VF and any subsequent 

19 medical intervention to def ibrillate, was thought 

2 0 merely to present unstructured electrical activity, and 
21 not to provide useful information. 

22 

24 from EKGs of the heart in VF focussed on the importance 

25 of the amplitude of the waveform defined using peak-to- 

26 trough differences in the EKG voltage, measured as 

27 either the greatest deflection occurring in a 

28 predefined time slot, or as the average peak- to- trough 

29 voltage measured over a given time interval. It has 

3 0 been shown that the VF amplitude is inversely related 

31 to time elapsed since collapse, is a crude predictor of 

32 defibrillation outcome, and is a better indicator of 
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1 survival to hospital discharge than the traditional 

2 variables described above . 
3 

4 However, recording the VF amplitude accurately is 

5 significantly problematical. The EKG voltage amplitude 

6 measured during VF is dependent on the direction of the 

7 main fibrillation vector and is influenced by a variety 

8 of factors including patient chest shape; electrode 

9 size; electrode location; and skin/electrode interface 

10 resistance- This number of variables makes this 

11 amplitude measure both unreliable and inaccurate. That 

12 is, although the amplitude of the waveform of an EKG 

13 recorded during VF is now recognised to be a crude 

14 predictor of the likely outcome of resuscitation of a 

15 patient in VF, it is not a reproducible marker of 

16 sensitivity to defibrillation, and lacks clinical 

17 usefulness. 
18 

19 In a further development, it is also known to use Fast- 

20 Fourier based transforms to generate a frequency 

21 spectrum of an EKG in VF to analyse the signal. The 

22 median frequency (MF) divides the area under the 

23 spectrum into two equal parts. Since this plot is \ 



24 derived from information in both the voltage and time 

2 5 domains, external variables such as lead placement have 

26 less effect on the results than the method of observing 

27 the amplitude. However, CPR produces artefacts in the 

28 recorded EKG signal and, since pausing CPR merely to 

29 obtain an EKG signal free of artefacts is likely to 

30 compromise resuscitation, these artefacts are 

31 necessarily included in this frequency measure, and 

32 detract from its usefulness. 



1 
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4 

1 

2 Thus the results of such signal analysis show some 

3 correlation with the likely outcome of resuscitation, 

4 but again lack sufficient sensitivity and specificity 

5 for clinical use. That is, this form of analysis has 

6 the disadvantage that, since the Fourier spectrum 

7 contains only globally averaged information, specific 

8 features in the signal are lost. 
9 

10 A method of accurate analysis of a surface EKG waveform 

11 recorded during VF would therefore be useful in 

12 understanding the pathophysiological processes in 

13 sudden cardiac death, and thus to produce a model for 

14 use: 
15 

16 in predicting the efficacy of therapy in individual 

17 cases; and 
18 

19 in determining the selection of the preferred course of 

2 0 primary, and alternative or adjunct therapies thus 

21 providing a means for individually tailored therapy for 

22 the specific patient needs 

23 " " 

24 to improve the success rate of resuscitation of 

25 patients presenting in VF. 
26 

27 Atrial fibrillation (AF) is a common cardiac arrhythmia 

28 in older people. Atrial fibrillation can be stopped by 

29 giving an electric shock to the patient under general 

3 0 anaesthetic (cardioversion) . However, many patient 

31 return to an AF rhythm soon after treatment. The 

32 technology detailed here may also provide a tool to 
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1 facilitate the clinical evaluation of AF exhibited in 

2 the electrocardiogram (EKG) so reducing the risk 

3 associated with general anaesthetic in patients where 

4 the applied therapy is likely to prove ineffective. 
5 

6 According to the present invention there is provided 

7 a method of decomposition of waveforms in a cardiac 

8 signal using wavelet transform analysis. 
9 

10 The method of the invention is non- invasive, accurate, 

11 and capable of delivering real-time information. 
12 

13 Preferably said method employs discretized wavelet 

14 transform analysis to process the EKG. 
15 

16 Preferably said method employs discretized continuous 

17 wavelet transform analysis to process the EKG. 
18 

19 Preferably said method comprises the steps of deriving 

20 the wavelet energy surfaces of an EKG signal; and 

21 plotting said wavelet energy surfaces against a 

22 location parameter b, and a scale parameter. The scale 

23 " parameter may be "dilation a"or" band pass frequency f h P n* 

24 

25 The method initially comprises the steps of connecting 

26 electrodes to the presenting patient; and sampling the 

27 analogue input signal to derive the cardiac signal . 
28 

29 Typically said method comprises the step of visually 

30 displaying the cardiac signal. 
31 
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1 Said method may display the distribution of energies 

2 within the cardiac signal. Said method may display 

3 coherent structures within the cardiac signal . 
4 

5 Said display may be by means of a contour plot. Said 

6 display may be by means of a surface plot. Preferably 

7 said method provides means to visualise the signal in 

8 real-time for clinical use. 

10 Preferably said method is applicable in the analysis of 

11 an EKG in ventricular fibrillation. 
12 

13 Said method may be applicable in the analysis of an EKG 

14 in ventricular fibrillation after the commencement of 

15 cardio-pulmonary resuscitation (CPR) . 
16 

17 The method may include the step of disassociating the 

18 component features of the temporal trace of a recorded 

19 EKG. Additionally or alternatively said method may 

20 include the step of temporal filtering of an EKG signal 

21 of a heart which is subject to CPR to disassociate the 

22 CPR signal from the heart signal . 

2 3 I "~~ 

24 Typically said method provides measurable 

25 characteristics for the estimation of the health of a 
2 6 heart in VF. Said method may provide measurable 

27 characteristics for the estimation of the health of a 

2 8 heart in AF. Said me may provide Typically said method 

29 provides measurable characteristics for the estimation 

30 of the health of a heart. 
31 
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7 

1 The method may provide measurable characteristics for 

2 the estimation of the time elapsed since the onset of a 

3 cardiac incident . 
4 

5 Typically said method provides measurable 

6 characteristics for the estimation of the health of a 

7 heart after commencement in CPR. 
8 

9" Said method may a "prediction 'for the outcome of 

10 a given therapeutic intervention and so aid the 

11 clinical decision making process. 
12 

13 Said method may provide a basis for individual, patient 

14 specific, protocols for therapeutic intervention. 
15 

16 The method may provide a guide to the optimal timing of 

17 defibrillation of a heart in VF. 
18 

19 Said method may include the step of constructing a 

2 0 damage index for reference purposes. Construction of 

21 said index might involve the development of a network 

22 classifier from a library of recorded data. Said 

2 3 n ^tw15r]^ cTa^^ a neural" Ti el iwork^ Said" ~ 

24 network classifier may comprise a wavelet network 

25 classifier. 
26 

2 7 Application of the method of the invention represents a 

2 8 significant advance in coronary care by providing a 

2 9 reliable predictor of the outcome of shocking a patient 

30 in VF. In addition, the development of an algorithm 

31 using the method of the invention gives the ability to 

32 predict shock outcome and to facilitate individual 
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1 patient therapy. The ability to provide patient 

2 specific therapeutic intervention is a priority in the 

3 advancement of currently applied medical protocols. 
4 

5 That is, as discussed above, in certain instances, 

6 after prolonged cardiac arrest preceding defibrillation 

7 pharmacological measures or CPR can increase the chance 

8 of successful resuscitation. Thus, employing the 

9 method to predict the outcome of shocking avoids futile 

10 defibrillation attempts which can even harm the heart, 

11 and can indicate the need for intervention, and 

12 influence the selection of the preferred type of 

13 intervention, to optimise the metabolic state of the 

14 heart prior to counter- shock . 
15 

16 The predictor algorithm developed using the method is 

17 being tested using a new generation of defibrillation 

18 devices that have the flexibility to allow easy 

19 prototyping of the new defibrillation algorithms. 
20 

21 According to a further aspect of the present invention 

22 there is provided a method of decomposition of 

23 waveforms in a cardiac signal "using matchin g pursuit 

24 algorithms . 
25 

2 6 According to a further aspect of the present invention 

2 7 there is provided an apparatus for decomposition of 

28 waveforms in a cardiac signal, said apparatus 

2 9 comprising wavelet transform analysis means. 
30 

31 Said apparatus may include means to display the 

32 distribution of energies within a waveform. 



_ PflgGB 00/ 0 1675 

9 

1 Said apparatus may include a monitor adapted to display 

2 decomposed waveforms. Said apparatus may be adapted 

3 for inclusion in an EKG apparatus. 
4 

5 According to a further aspect of the present invention 

6 there is provided defibrillation means adapted to 

7 operate in response to a signal generated by comparison 

8 of an EKG trace with decomposed waveform. 
9 

10 That is, the invention preferably provides a method of 

11 wavelet analysis of cardiac signals which provides 

12 structural information about the heart - whether the 

13 heart is healthy or not - and has significant 

14 advantages over fast Fourier transforms. 
15 

16 The invention may provide a display device in the form 

17 of a scrologram that provides real-time visualisation 

18 of a wavelet scalogram, showing the distribution of 

19 energies and coherent structures within the signal for 
2 0 use as guidance by a clinician. 

21 

22 The invention may further provide a data analysis tool, 

23 which" assists ~in ^ "That 

24 is, the derived data may indicate the optimum time to 

25 administer shock to the heart. The invention may 

26 provide a damage index, preferably in the form of an 

27 artificial neural network. 
28 

2 9 Preferably the invention provides dissociation of the 

3 0 component features of a temporal trace of a cardiac 

31 signal, which may for example be CPR, AF, or cardio- 

32 phonographic signals. 
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1 Embodiments of the invention will now be described by 

2 way of example only and with reference to the 

3 accompanying drawings in which: 



4 
5 

6 Figure la is a Mexican hat wavelet; 
7 

8 Figure lb is the real part of a complex Morlet 

9 wavelet; 
10 

11 Figure 2a is a schematic plot showing the dilation 

12 of a continuous wavelet; 
13 

14 Figure 2b is a schematic plot showing the 

15 translation of a continuous wavelet ; 
16 

17 Figures 3a to Figure 3e are the plots of the 

18 'investigation' of a sinusoidal signal by Mexican. 

19 hat wavelets of various sizes, showing the effect 
2 0 of translation of the wavelet along the signal 

21 (change in Jb) , and dilation of the wavelet (change 

22 in a) ; 

24 Figure 4a is the plot of five cycles of a sine 

2 5 wave of period P; 

26 

27 Figure 4b is the contour plot of T(a,b) against a 

2 8 and b for the sine wave of Figure 4a; 

29 

30 Figure 4c is the isometric surface plot of T(a,b) 

31 against a and Jb for the sine wave of Figure 4a; 
32 
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Figure 5a is the plot of a combination of two sine 
waves of period PI, and P2 , where PI = 5P2; 

Figure 5b is the contour plot of T(a.,b) against a 
and Jb for the sine wave of Figure 5a; 

Figure 5c is the isometric surface plot of T(a.,b) 
against a and b for the sine wave of Figure 5a; 

Figure 6a is an EKG trace of a pig heart in sinus 
rhythm; 

Figure 6b is a 2D energy scalogram associated with 
the EKG trace of Figure 6a; 

Figure 6c is a 3D energy scalogram associated with 
the EKG trace of Figure 6a; 

Figures 6d, 6e, 6f and 6g are the energy surface 
plots from four segments of an EKG signal 
subsequent to the onset of VF, showing the three 
dominant ridges A, B, and C appearing in the 
Itx_ans_f jorm^su th e 
onset of CPR after five minutes, associated with a 
gradual increase in passband frequency of the 
ridges A,B, and C; 

Figure 7a is an energy scalogram for a pig heart 
for the first seven minutes of ventricular 
fibrillation, indicating the initiation of CPR 
after five minutes; 
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1 Figure 7b is a schematic diagram of the salient 

2 features of the scalogram of Figure 7a; 
3 

4 Figure 7c is the smoothed plot of energy at the 

5 8Hz level in the scalogram of Figure 7a against 

6 time; 
7 

8 Figure 8a is a typical segment of an EKG trace of 

9 a pig heart in VF; 
10 

11 Figures 8b, 8c,. and 8d are the energy scalograms 

12 associated with the trace of Figure 8a; 
13 

14 Figure 9 is a screen shot of a real time viewer 

15 which shows the collected EKG data with its 

16 associated wavelet energy display in the form of 

17 its energy scalogram, where windows scroll to the 

18 right; 
19 

20 Figure 10a is a 7 second trace of human ECG 

21 showing a shock event; 
22 

"23 " ~ "F igure "1 Ob is a^scalogrim 'corresponding "to the ' ^^^^ 

24 trace of Figure 10a; 
25 

26 Figure 11a shows the proportion of energy in 

27 scalograms for 120 results (60 ROSC, and 60 

28 asystole) at 1.9 Hz after shocking ; 
29 

30 Figure lib shows the proportion of energy in 

31 scalograms for 120 results (60 ROSC, and 60 

32 asystole) at 9.3 Hz after shocking; 
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Figure 12a is a schematic representation of 
overlapping signal segments used in a neural 
network test study; 

Figure 12b shows the weights attributed by the 
Kohonen network to the 3 0 frequency levels used in 
the scalogram; 

Figure 13a is an aorta pressure trace ; 

Figure 13b shows the EKG for the same time period 
as the trace of Figure 13a; and 

Figure 13c is the scalogram associated with the 
trace of Figure 13a derived from the Morlet 
wavelet ; 

Figure 13d is a detail of the phase part of 
scalogram Figure 13c ; 



Figure 13e is the scalogram associated with the 
t rac e of "Figure "1" 3 a" ' der rved~f ~r ottT"t "Ke~Mex i~c aff "hat 



wavelet; and 

Figure 13 f demonstrates the correlation of aorta 
pressure pulse position with lines of zero phase; 

Figures 14a is the plot of an EKG trace. Figure 
14b is its associated phase at around 1.5Hz. 
Figure 14c is its energy scalogram. The 
correlation of zero phase at this lower frequency 
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1 and high frequency (low dilation) peaks is thus 

2 illustrated. 
3 

4 Figure 15a shows a 2 second segment of EKG taken 

5 from a patient with atrial fibrillation (AF) . 

6 Figure 15b shows the wavelet scalogram plot 

7 associated with this EKG. Figure 15c shows the 

8 corresponding modulus maxima of the scalogram of 

9 Figure 15b. 
10 

11 Figure 15d contains a 7 second segment of EKG 

12 exhibiting AF. Figure 15e is a trace of EKG 

13 temporal components with small amplitude. Figure 

14 15f shows the larger magnitude components i.e. the 

15 QRS and T waves. 
16 

17 Figure 15g is a plot of a two second x blow up' of 

18 part of the signal of Figure 15d; Figure 15h is a 

19 plot of a two second vblow up' of part of the 

20 signal of Figure 15e; and Figure 15i is a plot of 

21 a two second 'blow up' of part of the signal of 

22 Figure 15f. 

2 3 " ; 

24 Referring to the Figures, the present method employs 
2 5 the use of a wavelet transform to analyse a cardiac 

2 6 signal. 
27 

28 The method involves the decomposition of the signal. 

29 This decomposition is accomplished by utilising wavelet 

30 transforms to decompose the signal in wavelet space. 
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1 A key distinction between the Fourier analysis of an 

2 EKG signal and its analysis by means of a wavelet 

3 function is that, whilst the Fourier transform employs 

4 a sinusoid function, a wavelet function is localised in 

5 time. 
6 

7 The methodology for such decomposition may include 

8 discretized continuous wavelet transforms, orthonormal 

9 wavelet transforms of decimated construction, non- 
10 decimated wavelet transforms, wavelet packet transforms 
11 and matching pursuit algorithms. 

12 

13 Signal processing employing wavelet transform analysis 

14 allows simultaneous elucidation of both spectral and 

15 temporal information carried within a signal. Such 

16 processing can employ either continuous or discrete 

17 transforms. The choice of wavelet transform used for a 

18 particular signal processing application depends on 

19 factors such as speed of computation necessary, the 

20 shape of signal specific features, the frequency 

21 resolution required, and the statistical analysis to be 

22 performed . 

24 The preferred method employs the discretized continuous 

25 transform as it provides high resolution in wavelet 
2 6 space at lower frequencies. 

27 

28 This method thus employs the use of a discretized 

2 9 continuous wavelet transform to analyse a cardiac 

3 0 signal. 
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1 In particular, this method employs a wavelet transform 

2 as an interrogation tool for EKG signals of ventricular 

3 fibrillation. 
4 

5 A variety of wavelet functions are available, and the 

6 most appropriate is selected to analyse the signal to 

7 be investigated. 
8 

9 The wavelet transform of a continuous time signal, 

10 x(t) , is defined as: 
11 

12 T(a,b)=—^— f° x(t)m— -\lt equation 1 

w(a) Xm V a J 

13 

14 where g(t-b)/a) is the analysing wavelet function and 

15 1 ' denotes complex conjugate. w(a) is a scaling 

16 function usually of the form w(a)=a n where n is usually 

17 1 or 0.5, and x(t), in this application, is the single 

18 channel surface EKG time signal. The transform 

19 coefficients T(a,b) are found for both specific 

2 0 locations on the signal, b, and for specific wavelet 

2.1 dilations.,... a. T (^a,.jbA-„is_plo.t.ted_.against_„a.„and__Jb. .in 

22 either a surface or contour plot. 

23 

24 While other wavelet types may be employed the wavelets 

25 mainly used in this method are: the Mexican hat wavelet 

26 and the Morlet wavelet, examples of which are shown in 
2 7 Figure 1. 
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1 The wavelet can translate along the signal (change in 

2 i>) and dilate (change in a) . This is shown 

3 schematically in Figure 2 using a Mexican hat wavelet. 

4 Figure 3 illustrates the way in which a sinusoidal 

5 signal can be * investigated' at various locations by 

6 Mexican hat wavelets of various sizes. The numerical 

7 value of the convolution (equation 1) depends upon both 

8 the location and dilation of the wavelet with respect 

9 to the signal . 

10 Figure 3a shows a wavelet of similar 'size' to the 

11 sinusoidal waves superimposed on the signal at a Jb 

12 location which produces a reasonable matching of the 

13 wavelet and signal locally. From the Figure it is 

14 apparent that there is a high correlation between the 

15 signal and wavelet at this a scale and Jb location. 

16 Here, the cross correlation of the signal with the 

17 wavelet produces a large positive number T(a,b) . 

18 Figures 3b and 3c show details of the wavelet transform 

19 of a signal using a wavelet of approximately the same 

20 shape and size as the signal in the vicinity of Jb. 

21 -Figure 3b -shows a wavelet of similar -scale to the 

22 sinusoidal waveform located at maximum negative 

2 3 correlation. This produces a large negative T(a.,b) 

24 value. Figure 3c shows a wavelet of similar scale to 

25 the sinusoidal waveform located at a position on the 
2 6 time axis where near zero values of T(a,b) are 

27 realised. Figure 3d shows the effect on the transform 

28 of using the smaller a scale. It can be seen from the 

29 plot that the positive and negative parts of the 

30 wavelet are all in the vicinity of approximately the 
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1 same part of the signal, producing a value of T(a,b) 

2 near zero. Figure 3e shows that the same thing happens 

3 when using a much larger wavelet, since the wavelet 

4 transform now covers various positive and negative 

5 repeating parts of the signal, again producing a near 

6 zero value of T(a,b) . 

7 

8 Wavelet transforms are not usually computed at 

9 arbitrary dilations for isolated locations in the 

10 signal, but rather over a range of a and b. A plot of 

11 T(a, b) versus a and b for sinusoidal data using the 

12 Mexican hat wavelet is shown in Figure 4 . Two methods 

13 are then employed to plot T(a,b), namely a contour plot 

14 or scalogram as shown in Figure 4b, and a surface plot 

15 as shown in Figure 4c. At small and large values of a, 

16 the near zero values of T(a,b) are evident from the 

17 plots, but at values of a of the order of one quarter 

18 of the wavelength of the sinusoid large undulations in 

19 T(a,b) correlate with the sinusoidal forms of the 
2 0 signal . 

_2i _ — — 




22 Figure 5a shows two superpositioned sinusoidal 

23 waveforms, the first with period PI, the second with 

24 period P2 . PI = 5P2 . Figures 5b and 5c, the transform 
2 5 plots of the superimposed waveforms clearly show the 

2 6 two periodic waveforms in the signal at scales of one 

27 quarter of each period. Thus, Figure 5 clearly 

28 demonstrates the ability of the continuous wavelet 

29 transform to decompose the signal into its separate 
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1 frequency components. That is, this transform 

2 'unfolds' the signal to show its constituent waveforms. 

3 The contribution to the signal energy at a specific a 

4 scale and b location is proportional to the two- 

5 dimensional wavelet energy density function which is, 

6 in turn, proportional to the modulus of T(a,b) . 

7 

8 The method of the present invention thus involves the 

9 display of the transform as a contour plot. That is, 

10 the method is used to present information derived from 

11 an EKG trace of the heart in VF as a scalogram. The 

12 preferred form of presenting the information is as an 

13 energy scalogram, which presents the results as a plot 

14 showing the log of the wavelet energy coefficients, 

15 against the log of the bandpass centre frequency, fbpa 

16 of the wavelets for each time increment. The bandpass 

17 centre frequency is proportional to the reciprocal of 

18 the dilation value, a. This plot highlights small 

19 changes in amplitude over the scales of interest. The 
2 0 transform copes with repeating features in time with 

. 2 1 — s-hi-f-fei-ng— phase-? — making— i-t — appropri-a te— f or— real- t i me — — 

22 applications such as this. 
23 

24 That is, by performing continuous wavelet transform 

25 analysis on the ECG in VF, and then by producing an 

26 energy scalogram of the results, it is possible to 
2 7 unfold the signal in such a way that a previously 
2 8 hidden structure is apparent, in contrast to the 
29 apparently disorganised VF signal. 

30 
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1 The method then includes quantifying the wavelet 

2 decomposition. This wavelet decomposition provides 

3 both qualitative visual and measurable features of the 

4 EKG in wavelet space. 

5 

6 In practice, surface EKG tracings, recorded as soon as 

7 possible after the onset of VF, are analysed. 
8 

9 As a demonstration of the efficacy of the method, in an 

10 example of an experimental procedure utilising this 

11 method of analysis employing wavelet techniques, VF was 

12 induced in anaesthetised pigs via a pacemaker probe, 

13 using a 90V impulse at 60 Hz. All of the pigs remained 

14 in VF, untreated for a period of either 3 or 5 minutes. 

15 After this time, CPR commenced. The surface EKG 

16 (standard lead II) was recorded using needle 

17 electrodes. The EKG was sampled at 3 00 Hz using a 12- 

18 bit A to D converter. The method of the present 

19 invention was then performed using 32 EKG tracings 

20 recorded immediately after the onset of VF. 
21 

,22 F igure_ . 6 a„rjepres_e_nt_s___4 be ats of a pig he art i n _simis 



23 rhythm. Figures 6b and 6c shows the wavelet transform 

24 of the signal displayed in two and three dimensions 

25 respectively. 
26 

2 7 The QRS complex of the waveform is evident from the 

2 8 conical structures in Figure 6b converging to the high 

29 frequency components of the RS spike. The P and T 

30 waves are also labelled in the plot. The 3D landscape 

31 plot of Figure 6c shows the morphology of the signal in 
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1 wavelet space. In Figures 6b and 6c the continuous 

2 horizontal band (X) is associated with a frequency of 

3 1.7 Hz, the beat frequency of the sinus rhythm. The 

4 second band (Y) occurs at a frequency of approximately 

5 5.1 Hz, corresponding to the separation of the P-QRS-T 

6 components in time. At higher frequencies the P, QRS 

7 and T components are individually resolved according to 

8 their frequency makeup and temporal location. 
9 

10 Figures 6d to 6g show the energy surfaces for four 

11 segments of EKG signal subsequent to the onset of VF, 

12 namely: (6d) 0-60 s; (6e) 60-100 s; (6f) 210-240 s; 

13 and (6g) 260-360 s. 
14 

15 The morphology of the VF signal in wavelet space can be 

16 seen from the Figures to contain underlying features 

17 within a more complex surface topography. The most 

18 significant features are the dominant ridges that 

19 appear in the transform surface through time. 
20 

21 Figure 6f shows these ridges quite clearly. A high- 

2 2 energy ridge can be observed at around 10 Hz and two 

. 23 lc w^^nergy bands "can" be"observed at lower 

24 frequencies. These three ridges are labelled A, B and 

25 C, respectively, in the plot. Other ridges are also 
2 6 present within the scalogram. 

27 

2 8 The energy surface in Figure 6g contains the onset of 
29 CPR after 5 min of untreated VF . The institution of 

3 0 CPR is associated with a gradual increase in the 

31 passband frequencies of ridges A, B and C. This change 

32 in the composition of the VF signal reflects electrical 



RJT/GB 0 0 / 0 1 8 7 5 



22 

1 changes in the fibrillating myocardium associated with 

2 the onset of CPR. This is because CPR produces 

3 antegrade myocardial blood flow and thus improves the 

4 metabolic state of the tissues, temporarily reversing 

5 the otherwise progressive decline in high band pass 

6 frequency components of the EKG wavelet decomposition. 
7 

8 Figure 8a is a typical segment of an EKG trace of a pig 

9 heart in VF; Figures 8b, 8c, and 8d are the energy 

10 scalograms associated with the trace of Figure 8a. As 

11 clearly illustrated . by these diagrams the principle 

12 dilation (band pass centre frequency) component of the 

13 scalogram is approximately 10Hz. However, using said 

14 method it is also apparent that this component is not 

15 constant. It 'pulses' with a degree of regularity. This 

16 structure is previously unreported. 
17 

18 Figure 9 shows similar 'pulsing' in another porcine EKG 

19 signal. However, the structure is so pronounced that 
2 0 high energy, high frequency, intermittent components 

21 can be observed. These components have an occurrence 

22 frequency of the order of the original sinus rhythm: 

23 approximately 1 ,~7Hz ~ ~ 
24 

25 Figure 10a is a human EKG signal segment containing a 

26 shock event. Figure 10b is the corresponding wavelet 

27 scalogram. It is apparent from the scalogram of Figure 
2 8 10b that both high frequency spiking and an 

2 9 intermittent high-energy region are present in the 

30 vicinity of 10 Hz and also above 10Hz. 

31 
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1 The high frequency spiking is unique to the method of 

2 the present invention and is not visible using 

3 conventional Fourier techniques. The rich structure 

4 made visible within the EKG by the wavelet transform 

5 method is evident in the scalogram. 

6 It is clearly seen from the Figures that applying the 

7 wavelet transform to an EKG signal of VF demonstrates 

8 that this signal is a rich source of valuable 

9 information. That is, it produces a display showing 

10 real time visualisation of the distribution of energies 

11 and coherent structures within the signal for use by a 

12 clinician in the selection of treatment strategies. 

13 Using this method of analysis it is feasible to obtain 

14 real-time visual display of the EKG frequency 

15 characteristics in the wavelet domain during 

16 resuscitation. The scalogram produced provides 

17 information about the myocardium that is not available 

18 from a standard single channel surface EKG. 
19 

2 0 The wavelet scalogram decomposition can be displayed as 

21 a real time scrolling window, as shown in Figure 9. 

22 Thi s window is us e f ul ao an aid fo r clinic a l dec ision 

23 making. It can be used as a stand-alone tool, or as 

24 basis for on-line statistical analysis of the current 

25 state of a heart. 
26 

27 To produce the window, a MATLAB TM Rll application is 

28 used. Each EKG sample taken results in the updating of 

29 a FIFO (First In First Out) buffer, and the EKG plot of 

30 Figure 9a. The scalogram of Figure 9b is then shifted 
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1 to the right and clipped before the 'missing' new right 

2 hand data is calculated, using conventional matrix 

3 algebra, and filled. 
4 

5 This results in the two scrolling windows of Figure 9. 

6 The exponential ramp in the bottom right corner shows 

7 the compact support of the wavelet utilised at the 

8 given scale. 
9 

10 Higher resolution scalograms are achieved through 

11 implementation on higher specification machines, 

12 purpose built hardware, or application specific 

13 software with coding using a lower level programming 

14 language, such as C++. 

15 

16 CPR produces artefacts in the EKG signal. Additionally, 

17 this method delivers information the value of which is 

18 not degraded once the CPR artefacts are filtered from 

19 the EKG signal. 
20 

21 From examination of the scalograms shown in Figures 6g, 

-22 7a and 7b it can be seen that the VF signature and the — 

23 signature of the CPR artefacts occupy distinct areas of 

24 the scalogram, which permits their separation. 
25 

2 6 Known techniques such as the Modulus maxima method are 

27 now available to reduce the non-zero data points in the 

28 wavelet scalogram. This method reduces the topography 
2 9 of the scalogram surface to a series of ridges, thereby 



•CT/GB 0 0 / 0 1 6 7 5 



25 



1 considerably reducing the amount of data required to 

2 represent the signal in the wavelet space. 
3 

4 The modulus maxima obtained from a bandlimited signal 

5 with a wavelet of finite compact support in the 

6 frequency domain defines a complete and stable signal 

7 representation. 
8 

9 In this method/ temporal filtering of the original EKG 

10 signal to disassociate the CPR signature from the heart 

11 signal can either be done directly, using the wavelet 

12 energy scalograms , or indirectly through modulus maxima 

13 techniques. This allows the heart to be monitored 

14 without necessitating cessation of CPR to allow rhythm 

15 recognition. 
16 

17 Further to the above, the method may also be applied to 

18 patients suffering form atrial fibrillation (AF) as a 

19 means of disassociating the prevalent QRS and T waves 

20 from the remainder of the signal. 
21 

22 Wavelet decomposition of the ECG signal is performed 
"23" " using ah appropriate wavelet function . The modulus" 

24 maxima technique is used to encapsulate the scalogram 

25 information in a series of ridges. Filtering of the 
2 6 signal is then undertaken using the modulus maxima 

27 information and through reconstruction the clinically 

28 useful information is isolated from the signal . 
29 

30 Specifically, Figure 15a shows the wavelet transform 

31 decomposition of a 2 second segment of ECG taken from a 

32 patient with atrial fibrillation. Below the ECG trace 
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1 is a wavelet scalogram plot. The corresponding modulus 

2 maxima of the scalogram is plotted below the scalogram. 
3 

4 For example, Figure 15d contains a 7 second segment of 

5 ECG exhibiting AF. The signal has been partitioned 

6 using a modulus maxima ridge following algorithm. The 

7 modulus maxima ridges have been separated into large 

8 and small scale features by thresholding the signal at 

9 a predetermined wavelet scale. A blow up of part of the 

10 signal is given in the lower three plots in the figure: 

11 Figures 15g, 15h and 15i. The middle of these plots 

12 contains the partitioned signal with the QRS complex 

13 and T wave filtered out revealing regular, coherent 

14 features that appear at a frequency of approximately 

15 400 beats per minute, typical of AF. The lower plot 

16 contains the partition with the filtered out QRS and T 

17 waves. Although, a relatively simple modulus maxima 

18 technique was used in this pilot study whereby the 

19 modulus maxima lines were simple partitioned into two 

20 subsets, the ability of the technique to separate the 

21 signal into QRS and T waves and underlying atrial 

22 activity is evident from the results. It is known that 

23 the decay in ~ amplitude' of "a modijlug" ma xima ' 

24 corresponding to a signal feature can be a function of 

25 the scale of the wavelet. It is possible to use this 

26 property to separate the ridge coefficients into a 

27 noisy and coherent part. In this way, further 

28 differentiation of the modulus maxima information can 

29 be implemented within a more sophisticated algorithm. 

30 This will facilitate the further separation of 

31 background noise, QRS and T waves, and atrial activity. 
32 
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1 This method thus facilitates useful interpretation of 

2 previously unintelligible EKG signals. 

3- In patients presenting with uncoordinated rapid 

4 electric activity of the ventricle of heart, known as 

5 ventricular fibrillation (VF) , there is no effective 

6 pulse and myocardial blood flow ceases . Even the 

7 institution of optimal cardio-pulmonary resuscitation 

8 (CPR) of the patient does not achieve more than 3 0% of 

9 the normal cardiac output. Ischaemia during cardiac 

10 arrest leads to a rapid depletion of myocardial high- 

11 energy phosphates, deterioration of transmembrane 

12 potentials, and disruption of intracellular calcium 

13 balance. Paradoxically, the myocardium in VF has 

14 supranormal metabolic demands. For this reason 

15 resuscitation attempts become less likely to succeed 

16 with the passage of time, and electrical def ibrillating 

17 shocks increasingly result in asystole or EMD. 
18 

19 After prolonged cardiac arrest, the use of 

20 pharmacological measures or CPR before attempting 

21 defibrillation may increase the chances of successful 

2 2 res usc itati on. Th is invention provide s a robu st and 

23 reliable method of analysis of the state of the 

24 myocardium in VF that prevents attempts to defibrillate 

25 at times that are unlikely to be successful, or even 

26 harmful to the heart. This method also provides an 

27 indication of the best way in which to optimise the 

28 metabolic state of the heart prior to counter- shock . 
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1 The method includes steps to establish a standard 

2 against which to evaluate collected data in a 

3 particular incidence. 
4 

5 The method further employs use of measurable signal 

6 characteristics derived from the position and amplitude 

7 of features in the scalogram to estimate both the 

8 condition of the myocardium, and downtime of the 

9 subject while in VF. 
10 

11 The method thus provides for optimal treatment of the 

12 heart in VF, so fulfilling specific patient needs, by 

13 therapeutic intervention, if appropriate. 
14 

15 An energy scalogram such as that shown in Figure 7 

16 displays three distinct bands, labelled A, B, C. It is 

17 possible to derive quantifiable measures using 

18 correlations between the location and energy content of 

19 the bands . 
20 

21 Band A of Figure 7b represents the dominant energy band 

22 seen in the scalogram of Figure 7a, and corresponds to 
23~ the~tacHyca"rdi'c~b^atri"hg of VF7 However" the scaTo'gram 
24 is much more informative in that it also shows, as 

2 5 bands B and C, the behaviour of other frequency 

2 6 components of the signal which were previously 

2 7 unreported . 
28 

2 9 Figure 7a shows a 2D energy scalogram. It includes the 

30 first 5 minute period of VF, followed by a 2.5 minute 

31 period of CPR. The onset of CPR is clearly identified 

32 by the distinct horizontal dark band in the lower right 
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1 quadrant of the Figure. Over the first 5 minute 

2 period, three bands, labelled A, B # C, can be clearly 

3 seen in the scalograms. These bands correspond to the 

4 ridges of Figures 6d to g. The increase in the 

5 frequency components of these three bands after the 

6 onset of CPR is evident in the plot. Bands B and C 

7 follow trajectories similar to each other in the 

8 scalogram, reducing in frequency over time. Band A, 

9 however, moves independently of the other two. 

10 Initially Band A increases, then it decreases to a 

11 local minimum value at approximately 70s. Between 70 

12 and 160s it increases relative to Bands B and C. 

13 Finally, it decreases until the start of CPR after 

14 300s. The same pattern was present in all 32 pig EKG 

15 traces of the experiment . 
16 

17 Obvious increases in the passband frequency of all 

18 three bands are observed in the scalogrram after the 

19 onset of CPR. For some of the signals studied this 

20 increase in band C is masked by the dominant CPR band, 

21 and thus cannot be seen in the scalogram. 
22 

23 " Fiqure~~7b~provi~de~s~^ of "the "salient " ' 

24 features contained within the scalogram plots, where to 

25 is immediately after the onset of VF; t2 is the start 
2 6 of CPR; and t3 is the end of the analysis. Figure 7c 

27 shows the relative proportion of energy contained in 

28 the scalogram in the 5 to 12 Hz region through time. 
2 9 There is an obvious decay in the relative energy 

30 associated with this region which is associated with 

31 the breakdown of co-ordinated activity in the heart. 
32 
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1 The steps of the method of the present invention 

2 described above establish that during the course of VF 

3 there is a reduction in the proportion of energy within 

4 the dominant frequency band indicated in Figure 7c . 

5 This dominant frequency band, Band A in Figure 7a, is 

6 demonstrated to be approximately 10 Hz for pig VF. 
7 

8 The energy within this band changes rapidly. This is 

9 illustrated by the 'pulses' in Figures 8,9,10. 
10 

11 The Figures 6,7,8,9,10 show that applying the wavelet 

12 transform to an EKG signal of VF demonstrates that this 

13 signal is a rich source of valuable information. 
14 

15 The underlying hypothesis of the method of the present 

16 invention is that the scalogram associated with an EKG 

17 correlates to the state of the myocardium as it decays 

18 subsequent to the onset of VF. 
19 

2 0 The method uses the information contained in the energy 

21 scalogram associated with an EKG to predict the likely 

22 success of clinical intervention, namely shocking. 

24 It is therefore possible to develop a wavelet transform 

2 5 based tool for the prediction of shock outcome during 

26 ventricular fibrillation by: 
27 

28 1. collecting and collating data from sets of 

2 9 archived EKGs recorded from humans in VF where 

30 attempts to resuscitate by shocking were made; and 

31 

32 2. developing a classifier for reference purposes. 
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1 

2 Figure 11 is a classification of the shock outcome in 

3 either asystole or a rhythmic response using a 

4 relatively simple statistical analysis. The experiment 

5 yielding the results to compile these Figures involved 

6 use of the lead II outputs of standard three lead EKGs 

7 of 120 patients in VF . Each trace is of three second 

8 duration sampled at 100 Hz. Of these patients, 60 

9 returned to sinus rhythm while the other 60 
10 deteriorated to asystole, post shock. 

11 

12 Each trace was decomposed into an associated wavelet 

13 transform from which its energy scalogram was 

14 generated. The volume under this surface was then 

15 normalised to render the results independent of signal. 

16 amplitude, but instead the result of the relative 

17 wavelet constituents of the signals. The log of the 

18 mean values at each dilation (band centre frequency) 

19 for each was then recorded. Figures 11a and lib show 

20 the distribution of energies in a lower frequency band 

21 (1.9 Hz) and at the 9.3 Hz band. Clearly, through 

22 visual inspection, it is apparent that the proportion 

23 of- -energies -around- -the— 1-0- -Hz— band- is- higher -for 

24 successful defibrillation attempts. 
25 

2 6 The method then extends to apply neural techniques to 

27 analysis of wavelet pre-processed EKG signals. 

28 

29 A pilot study conducted to determine the feasibility of 

30 using artificial neural techniques to provide a tool to 

31 predict the outcome of defibrillation during VF used 

32 eight human EKG trace segments containing shock events. 
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1 In these cases, the result of shocking was unequivocal 

2 - four patients returned to VF, and four experienced 

3 return of spontaneous circulation (ROSC) . 
4 

5 The traces were transformed using the Morlet wavelet, 

6 and energy scalograms containing thirty frequency 

7 levels were produced. This was then split into eight 

8 overlapping sections as shown in Figure 12a, each of 

9 2 00 points (2/3 seconds duration) . These 200 location 

10 points were subsampled down to 50 to give eight 

11 scalograms for each trace of 50 x 30 elements. The 

12 volume under the energy scalograms were normalised and 

13 the patterns fed into a 'winner take all' Kohonen 

14 network with two output units and built in conscience 

15 (to avoid local minima) . That is, the network was 

16 asked to group the 64 input patterns into two classes. 

17 All but ten outputs were collectively classified 

18 correctly giving a mean pattern error of 0.156 (against 

19 0.5 average pattern error expected from random inputs) . 

20 _ 

21 Since this is a vector quantisation method (VQM) it was 

22 possible to identify how the network differentiates the 
2 3 pa t: terns thr ouqtrTnspe c trion~o f ~i~t~s~colin~e^ 

24 The weights from each location position across all 

2 5 scales in the network are approximately the same, which 

2 6 means that there are no markers with which to 

27 synchronise the different pre-processed traces. This 

2 8 confirms that this neural network is too simple for 
29 this purpose. That is the network is not equipped to 

3 0 'consider' the relative phase of each input pattern. 
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1 Figure 12b shows the weights for the 'success' (ROSC) 

2 and x failure' (VF) to the output units from the first 

3 two time slices across all scales. The weights 

4 indicate the classes are differentiated by the 

5 proportion of energy in the lower scales, which can be 

6 seen when compared with Figure 11. 
7 

8 Although the above described method indicates the 

9 slight drop in the dominant frequency expected, the 

10 drop is very marginal which leads to the conclusion of 

11 the lack of competence of previously proposed methods 

12 as a defibrillation success predictor. 
13 

14 In summary, a library of human ECG data containing data 

15 sets of human VF with attempts to resuscitate by 

16 shocking is used as a database. This database is 

17 extended to include data sets containing various 

18 methods of shocking including, for example, biphasic 

19 shocking. The biphasic shock waveform has resulted in 

20 an increased proportion of successful defibrillation 

21 attempts and is set to become the standard treatment 

22 for cases of VF. 
23 

24 In one example, the recognised outcomes are defined by 

2 5 trace components of the post -shock window lasting until 

2 6 next shock (if present) . If the ratio of the given 

27 rhythm exceeds 10% of the total window length the 

28 rhythms are prioritised according to the sequence: 
29 

30 Class Rhythm Ratio 

31 

32 1 Pulse (SVR) +10% 
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12 No pulse (EMD) +10% 

2 3 Isoelectric (Asystole) +10% 

3 4 VF +10% 
4 

5 

6 Class 5 is the class of VF preceding shocks where VF 

7 re-establishes itself within 5 seconds following the 

8 shock (i.e. no change) . The VF in all the other 

9 classes were non-VF in this period. 
10 

11 Wavelet analysis of this information in accordance with 

12 the method of the invention is then performed to: 
13 

14 construct a wavelet visualisation of the signal - 

15 usually by plotting wavelet energy surfaces against the 

16 location parameter b and the inverse of the dilation 

17 parameter a; 
18 

19 provide measurable characteristics of the signal for 

20 estimation of downtime of the patient; 
21 

22 provide measurable characteristics of the signal for 

23 determining"the "health of "the""heart~'post" CPR;"and — — 

24 

2 5 to construct energy scalogram devised for the method - 

2 6 which uses the energy density function and the 

27 reciprocal of the wavelet a scale for use as a 

28 predictor tool. 
29 

30 As described above it is possible to use artificial 

31 neural network based techniques to develop such an 

32 indication of the state of myocardium. In the 
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1 alternative, it is possible to classify the wavelet 

2 scalogram through multilayered feedforward network 

3 types . 
4 

5 The method may include the development of a modulus 

6 maxima algorithm tool for the preprocessing of ECG 

7 prior to its input into a neural network classifier. 
8 

9 Using this technique improves network performance 

10 whether this data is further encoded, or presented as a 

11 whole, larger, sparse matrix as a pattern in the input 

12 space. 
13 

14 This method therefore utilises the generalisation 

15 properties of a feed forward multi-layer network to 

16 predict the likelihood of defibrillation success from 

17 the wavelet transform of the EKG traces. This multi- 

18 layer network with its relatively simple dynamics, when 

19 combined with wavelet pre-processing, has proved itself 

20 a useful tool as a universal approximator. 
21 

22 The classes of multi-layer network types of use in this 

2 3 method—are-: 

24 

2 5 • Mult i- layered feed forward (MLFF) neural networks 

26 with back propagation training and monotonic 

27 activation functions; and 

2 8 • Radial Basis Neural Networks (RBNN) as have 

29 previously been successfully applied to the denoising 

30 of medical Doppler ultrasound signals with wavelet 

31 preprocessing. 
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1 

2 As described above, the method involves the 

3 decomposition of EKG signals into a complete basis set 

4 defined by the wavelet shape and other parameters by 

5 salient basis functions of a different basis set, 

6 converged upon through regression techniques (sigmoid 

7 in the case of multilayer neural networks, Radial basis 

8 etc) . 
9 

10 These regression techniques can also be used to 

11 construct a wavelet basis function set directly. 
12 

13 Methodologies for restricting the search space of the 

14 wavelet basis functions considered are known. Whilst 

15 this wavelet network has been shown to be effective for 

16 chaotic time series prediction, its implementation 

17 involves the use of wavelet frames of a decimated, 

18 dyadic, construction. The method of the present 

19 invention may employ continuous wavelet networks 

20 spanning a redundant wavelet basis which, although 

21 computationally more expensive, overcomes the time 

22 invariance constraint and the limited size of input 
-23 space~assocrated"wrth-use"of"wave-l-et— f rames . 

24 

2 5 The method may use conventional gradient decent methods 

26 to produce a single layer wavelet classifier. 

27 

2 8 These wavelet networks may be further employed as part 

2 9 of a multilayer system as a non-parameterised estimate 

3 0 of the original trace for input to further hidden 
31 layers. 
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1 The network type of choice for the automated prediction 

2 system of the method is selected on the basis of its 

3 sensitivity and selectivity in correctly classifying 

4 successful defibrillation outcomes in test set data, 

5 since this is most clinically useful. 
6 

7 Thus experimental comparison of the three techniques 

8 demonstrates the efficacy of the wavelet transform 

9 technique . 
10 

11 The nature of underlying atrial activity can also be 

12 determined from wavelet decomposition of the EKG 

13 signal. The wavelet function gives information 

14 regarding the amplitude and, where appropriate, phase 

15 of the transformed signal. It is known that pressure 

16 readings taken from the aorta correlate to forms of 

17 atrial activity within the heart. Areas of localised 

18 high energy contained within the scalogram can be 

19 demonstrated to correlate with these pressure readings. 

20 This experimental result is extrapolated to mean that 

21 areas of localised high energy contained within the 

22 scalogram correlate with forms of atrial activity 

-23 wi thin— the—hear t . — 

24 

2 5 Figure 13a shows the aorta pressure, Figure 13b the EKG 

26 trace, for the same time period as Figure 13a, and 

27 Figure 13c shows the scalogram for the EKG of Figure 

28 13b. It is apparent that there is an increase in 

29 energy in the system during an atrial pulse, indicated 

30 by the dark blotches occurring in the scalogram at an 

31 f bpc of around 10 Hz. There is a frequency component 

32 between 1 and 2 Hz . As shown in Figure 13d, which 
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1 highlights the phase of the scalogram between 1 and 2 

2 Hz, it is apparent hat the lines of zero phase are in 

3 alignment with the atrial pulse. 
4 

5 In a further scalogram, shown in Figure 13e, produced 

6 by using the Mexican hat wavelet transform which is 

7 real and has better temporal resolution, but worse 

8 frequency resolution than the complex scalogram of 

9 Figure 13c, it is demonstrated that positive high 

10 amplitude components are shown at the same positions 

11 for scales of between 1 and 2 Hz, thus reinforcing the 

12 findings extrapolated from Figure 13c. That is as 

13 shown in Figure 13f, the lines of zero phase correlate 

14 with the pulse position. 
15 

16 The lines of zero phase within the 1 . 8Hz frequency band 

17 also align with regular peaks in the scalograms, as 

18 shown in Figures 14a, 14b & 14c. This links the 

19 presence of the 1.8 Hz band with the observed peaks at 
2 0 higher frequencies. This correlation between the 1.8 

21 Hz band and the aorta pressure pulse suggests atrial 

22 activity is present. 

___ 

24 In a further application of the method, means for 

25 identifying the optimum timing for application of the 

26 defibrillation shock can be extrapolated from the 

27 pulsing identified by the wavelet technique and shown 

28 in Figures 8, 9, 10, and 14, by comparison with traces 

29 of attempts at defibrillation which initially fail but 

30 are subsequently successful. 
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1 

2 Thus, any data sets, in the above, that correspond to 

3 multiple shocking of the same patient, where 

4 defibrillation has been repeatedly attempted are 

5 considered separately since these traces hold important 

6 information. 
7 

8 The pilot study detailed above used Morlet wavelet 

9 based energy scalogram decomposition of signal segments 

10 immediately prior to shocking. A full parametric 

11 wavelet study of the method determines the optimum 

12 method. 
13 

14 The method includes the development of a classifier 

15 using the wavelet transform analysis. 
16 

17 Various types of neural network classifier are 

18 achievable using this method. 

19 

20 The linkage of shock timing to the phase information of 

21 wavelet components allows for increased defibrillation 

22 success and reduced shock energies. The wavelet- 

"23 derived information can also be employed to predict the 

24 likelihood of shock success, preventing futile or 

25 harmful defibrillation attempts, and providing a 

26 predictor of an optimal resuscitation strategy or 

27 strategies. 
28 

2 9 This method demonstrates the utility of the wavelet 

30 transform as a new method of EKG signal analysis during 

31 VF. It provides a robust, real-time solution to the 
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1 problem of useful monitoring of the myocardium during 

2 resuscitation. 
3 

4 When compared with conventional statistical methods, 

5 such as fast Fourier transforms, it is seen that the 

6 temporal resolution of the wavelet technique gives a 

7 scalogram which better describes the non- stationary , 

8 intermittent, nature of the EKG trace to be analysed, 

9 and gives a method of greater predictive effectiveness 

10 than is already known. The effectiveness criteria for 

11 the networks of the method of the present invention are 

12 based upon their sensitivity and selectivity in 

13 correctly classifying successful defibrillation 

14 outcomes from test data sets. 
15 

16 Although this description refers to wavelet transform 

17 analysis, this term is to be construed to include 

18 matching pursuit algorithms and similar analysis 

19 techniques. 
20 

21 Modifications and improvements can be made to the above 

22 without departing from the scope of the invention. 
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1 CLAIMS 
2 

3 1. A method of decomposition of waveforms in a 

4 cardiac signal using wavelet transform analysis. 
5 

6 2 . A method as claimed in Claim 1 comprising the step 

7 of employing discretized wavelet transform 

8 analysis to process the said waveform. 
9 

10 3. A method as claimed in Claim 1 comprising the step 

11 of employing discretized continuous wavelet 

12 transform analysis to process the cardiac 

13 waveform. 
14 

15 4 . A method as claimed in any preceding claim 

16 comprising the steps of deriving the wavelet 

17 energy surfaces of an electrocardiogram (EKG) 

18 signal; and plotting said wavelet energy surfaces 

19 against a location parameter b, and a scale 
2 0 parameter. 

21 

22 5. A method as claimed in Claim 4 wherein said scale 

— 23 parameter-is— di-l-atron-~aT~ — 

24 

25 6. A method as. claimed in Claim 4 wherein said scale 

2 6 parameter is band pass frequency fb pc - 

27 

28 7. A method as claimed in any preceding claim 

29 comprising the initial steps of connecting 

30 electrodes to a presenting patient; and sampling 

31 the analogue input signals recorded to derive the 

32 cardiac signal. 
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1 8. A method as claimed in any preceding claim 

2 including visually displaying the cardiac signal. 
3 

4 9. A method as claimed in any preceding claim 

5 including visually displaying the distribution of 

6 energies within the cardiac signal . 
7 

8 10. A method as claimed in any preceding claim 

9 including visually displaying coherent structures 
10 within the cardiac signal. 

11 

12 11. A method as claimed in any preceding claim 

13 including visually displaying the signal in real- 

14 time for clinical use. 
15 

16 12. A method as claimed in any preceding claim 

17 comprising the step of constructing a contour plot 

18 to display the decomposed waveform obtained. 
19 

20 13. A method as claimed in any preceding claim 

21 comprising the step of constructing a surface plot 

22 to display the decomposed waveform obtained. 

"23 — — 

24 14 . A method as claimed in any preceding claim 

25 comprising the step of constructing a 2D or a 3D 

2 6 energy scalogram to display the decomposed 
27 waveform obtained. 

28 

29 15. A method as claimed in any preceding claim 

3 0 including the step of disassociating the component 
31 features of the temporal trace of a recorded EKG. 
32 
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1 16. A method for the analysis of an EKG of a heart in 

2 ventricular fibrillation including the method as 

3 claimed in any preceding claim. 
4 

5 17. A method for the analysis of an EKG of a heart in 

6 ventricular fibrillation after the commencement of 

7 cardio -pulmonary resuscitation (CPR) including* the 

8 method as claimed in any of Claims 1 to 15. 
9 

10 18. A method as claimed in Claim 17 including the step 

11 of temporal filtering of the EKG signal of a heart 

12 that is subject to CPR to disassociate the CPR 

13 signal from the heart signal . 
14 

15 19. A method as claimed in Claim 17 or Claim 18 using 

16 wavelet energy scalograms . 
17 

18 20. A method as claimed in Claim 17 or Claim 18 using 

19 ridge following techniques 
20 

21 21. A method as claimed in Claim 20 wherein said ridge 

22 following techniques are modulus maxima 

23 t'echni"qu'e"s" 
24 

25 22. A method for the estimation of the health of a 

26 heart in VF including the method of any of Claims 

27 1 to 15 to provide measurable characteristics. 
28 

29 23. A method as claimed in Claim 22 wherein said 

30 measurable characteristics are used to provide an 

31 estimate of the time elapsed since the onset of a 

32 cardiac incident. 
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1 24. A method as claimed in Claim 22 wherein said 

2 measurable characteristics are used to provide an 

3 estimate of the health of a heart after 

4 commencement of CPR. 
5 

6 25. A method as claimed in any of Claims 22 to 24 

7 wherein said measurable characteristics are used 

8 to predict the outcome of a given therapeutic 

9 intervention . 
10 

11 26. A method as claimed in any of Claims 22 to 25 

12 wherein said measurable characteristics are used 

13 to provide a guide for the optimal timing of 

14 defibrillation of a heart in VF . 
15 

16 27. A method for the analysis of an EKG of a heart in 

17 atrial fibrillation including the method as 

18 claimed in any of Claims 1 to 14 . 
19 

20 28. A method as claimed in Claim 27 including the step 

21 of partitioning the signal to provide separate 

22 traces of QRS and T waves, and/or atrial activity 

23 andy or ~bac kg round~noise~ 
24 

25 29. A method as claimed in any preceding claim 

2 6 including the step of constructing a damage index 

27 for reference purposes. 

28 

29 30. A method as claimed in Claim 29 wherein 

30 construction of said index includes the step of 

31 developing network classifier from a library of 

32 recorded data. 
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1 31. A method as claimed in Claim 30 wherein said 

2 network classifier developed is a neural network. 
3 

4 32. A method as claimed in any of Claims 29 to 31 

5 wherein said network classifier developed is a 

6 wavelet network classifier. 
7 

8 33. A method of decomposition of cardiac waveforms 

9 using matching pursuit algorithms. 
10 

11 34. Apparatus for decomposition of waveforms in a 

12 cardiac signal, said apparatus comprising wavelet 

13 transform analysis means. 
14 

15 35. Apparatus as claimed in Claim 34 including means 

16 to display the distribution of energies within a 

17 waveform. 
18 

19 36. Apparatus as claimed in Claim 34 or Claim 35 

2 0 including a monitor adapted to display decomposed 

21 waveforms. 

22 

23 3 7 . Apparatus^ a's~craimed^any"of~Cl"aims~3"4~to~~3 6~ 

24 adapted for inclusion in an EKG apparatus. 
25 

26 38. Defibrillation means adapted to operate in 
2 7 response to a signal generated by comparison of an 

2 8 EKG trace with decomposed waveform obtained by the 

29 method of any of Claims 1 to 33. 
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1 39. A method as described in any of Claims 1 to 33 

2 with reference to or as shown in the accompanying 

3 drawings . 
4 

5 40. Apparatus as described in any of Claims 34 to 38 

6 with reference to or as shown in the accompanying 

7 drawings . 
8 

9 
10 
11 
12 
13 
14 
15 
16 
17 
18 
19 
20 
21 
22 
23 
24 
25 
26 
27 
28 
29 
30 
31 
32 
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1 ABSTRACT 
2 

3 A method of analysis of medical signals which uses 

4 wavelet transform analysis to decompose cardiac 

5 signals. Apparatus for carrying out the method, and 

6 cardiac apparatus adapted to employ the method are also 

7 described. 
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